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We report new results on identified (anti)proton and charged pion spectra at large transverse 
momenta (3<pt<10 GeV/c) from Cu+Cu collisions at ^/snn = 200 GeV using the STAR detector 
at the Relativistic Heavy Ion Collider (RHIC). This study explores the system size dependence of 
two novel features observed at RHIC with heavy ions: the hadron suppression at high-pr and the 
anomalous baryon to meson enhancement at intermediate transverse momenta. Both phenomena 
could be attributed to the creation of a new form of QCD matter. The results presented here bridge 
the system size gap between the available pp and Au+Au data, and allow the detailed exploration for 
the on-set of the novel features. Comparative analysis of all available 200 GeV data indicates that 
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the system size is a major factor determining both the magnitude of the hadron spectra suppression 
at large transverse momenta and the relative baryon to meson enhancement. 



PACS numbers: 25.75.-q 

Differential studies of identified particle production in 
nucleus-nucleus collisions provide an experimental means 
to probe the different stages of the collision evolution and 
explore the properties of the created medium. Spectral 
measurements at high transverse momenta are of special 
interest, for the following reasons. In elementary colli- 
sions, hard partonic scatterings are known to produce jets 
of particles originating from the fragmentation of a high- 
Pt quark or gluon. The spectral distributions of particles 
in transverse momentum from such interactions are mea- 
sured experimentally and are reasonably well understood 
in terms of Next-to-Leading Order (NLO) pQCD calcu- 
lations [l[ . These hard scatterings occur in heavy- ion col- 
lisions as well, but their resulting distributions are found 
to be modified due to interactions with the medium and 
the resulting energy loss. Thus, understanding modifica- 
tions to the high-px particle distributions is an impor- 
tant step towards understanding the partonic energy loss 
mechanisms within the medium (2j. 

To study the effects of parton-medium interaction on 
particle production in heavy-ion collisions we compare 
the production cross-sections measured in AA to the 
equivalent measurements in pp collisions. Following the 
expectation that the particle production in heavy-ion col- 
lisions at high-px is determined by the number of binary 
nucleon-nuclcon inelastic collisions we define the nuclear 
modification factor, Raa, as the ratio (Eq. (Q}) of parti- 
cle yields measured in AA to the cross-sections measured 
in pp collisions scaled by the corresponding number of 
independent nucleon-nucleon collisions N^. We obtain 
^bin f rom a Monte Carlo Glauber model calculation Q . 
For the unmodified particle production in AA collisions 
Raa is exactly unity, whilst Raa<1 indicates suppres- 
sion and Raa > 1 enhancement. 

cr%% d^NAA/dydpT 

N££ d 2 ° pp /dydp T ■ W 

The Raa measured in d+Au and peripheral Au+Au col- 
lisions exhibits an enhanced particle production which 
is believed to occur due to multiple nucleon scatterings 
within the colliding nuclei. This "initial" state effect 
is known as the Cronin effect [3, 0, Meanwhile, in 
central Au+Au collisions, Raa at high-pr indicates that 
the particle production is strongly suppressed (by about 
a factor of 5) [H, ||. This "final" state effect has been 
attributed to the partonic energy loss in an opaque col- 
ored medium Q. However, neither of the two effects is 
sufficiently understood and both require further exper- 
imental and theoretical study. The differential analysis 
presented here explores the system size effects on parton 
propagation through the medium to further evaluate the 



mechanisms of parton/mcdium interactions. 

To provide additional constraints and systematic un- 
derstanding of the measurements in very light (d+Au) 
and heavy (Au+Au) collision systems we present the key 
studies at the intermediate (Cu+Cu) system at the same 
incident energy (^snn = 200 GeV), bridging the gap be- 
tween the two extremes. These measurements may pro- 
vide quantitative understanding of the partonic energy 
loss and its system size dependence. In addition, it is ex- 
pected that the identified particle measurements provide 
information on color-charge effects within the mechanism 
of jet quenching. Although experimental discrimination 
between quark and gluon jet fragmentation on event-by- 
event basis is difficult, it can be addressed statistically by 
exclusive analysis of proton (or baryon) and pion (meson) 
production. We are utilizing the idea that the baryon to 
meson ratio is found higher in gluon jets compared to 
quark jets [I]- Identified proton and pion measurements 
from pp collisions concur with this picture [H,Q, as well as 
direct measurements of baryon and meson production in 
quark and gluon jets [To| . Thus, identified particle mea- 
surements at high-pT can then be used to analyze gluon 
and quark propagation through the medium and to pro be 
the color-charge differences of energy loss 0, H, ITT1. Il2j . 

Additionally, systematic studies of identified parti- 
cle production in Cu+Cu can shed new light on the 
anomalous enhancement of baryons with respect to 
mesons observed at intermediate transverse momenta 
(2<pt<6 GeV/c) in Au+Au collisions. This enhance- 
ment is not consistent with the extrapolated values from 
the measurements in pp collisions and cannot be ex- 
plained by cold nuclear matter effects. At present, the 
preferred baryon over meson production at intermediate 
Pt could be described by two very different considera- 
tions. The first model assumes coalescence and recombi- 
nation, which demands a shift of baryon yields to higher 
momenta relative to meson yields [l3j . The second model 
evokes an interplay of the flow effects in the radially ex- 
panding medium with the jet fragmentation [Til ] . 

In this paper identified (anti)proton and charged pion 
spectra are systematically explored with regard to the 
system size with the smaller Cu+Cu colliding system at 
y/sNN= 200 GeV. The centrality dependence of high- 
Pt hadron production and the pt dependence of baryon 
to meson ratios in Cu+Cu data are compared to the 
Au+Au system as well as to pp collisions at the same en- 
ergy. This allows gaining a greater understanding of pe- 
ripheral collisions. The size of Cu nuclei is ideally suited 
to explore the turn-on of the high-pT suppression bridg- 
ing the gap between pp, d+Au and peripheral Au+Au 
data in terms of system size and nuclear matter. 

The Cu+Cu data used in this analysis were recorded 
by the STAR experiment during Run 5 at RHIC. Here, 
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FIG. 1: Transverse momentum spectra of pions (a) and protons (b) produced in Cu+Cu collisions at y / sjvjv=200 GeV. Data 
are presented for four centrality classes: 0-10%, 10-20%, 20-40% and 40-60%. Closed and open symbols are used for particles 
and antiparticles, respectively. For clarity, data are separated by powers of four. 



the minimum bias trigger was based on the combined sig- 
nals from the Beam-Beam Counters at forward rapidity 
(3.3<|?7|<5.0) and the Zero-Degree Calorimeters, located 
at ±18 m from the nominal interaction point [la ]. In to- 
tal, 23 million events comprise this data set. Based on 
the charged track multiplicity recorded in the Time Pro- 
jection Chamber (TPC) and Glauber MC model calcula- 
tions, the data are divided into four centrality bins corre- 
sponding to 0-10%, 10-20%, 20-40% and 40-60% of frac- 
tional cross-section (a/a geom ) bins. In order to remove as 
many background tracks as possible, tracks which inter- 
cept the measured collision vertex within 1 cm (distance 
of closest approach) were retained with a minimum of 25 
(out of 45) TPC trajectory points forming each track. 

Within the STAR experiment, particle identification 
at Iow-pt is attained by use of the ionization energy loss 
(dE/dx) in the TPC [lg]. For low momentum particles, 
below 1 GeV/c, a clear mass separation is observed al- 
lowing the identification of tt , K and (anti)protons. In 
the intermediate-pT region (l<pr<3 GeV/c) the TPC is 
no longer directly usable by itself, as all particles, inde- 
pendent of mass, are minimum ionizing. For the pur- 
pose of this paper we identify pions, kaons, protons and 
anti- protons at higher momenta (pr> 3 GeV/c) on a 
statistical basis utilizing the relativistic rise of the ion- 
ization energy loss in the TPC. For a given slice in trans- 
verse momentum, a distinctly non-singlc-Gaussian shape 



is observed, discussed in detail in [l7| . representing the 
normalized deviations from different energy loss trends 
of tt, K and protons. The quantity used to express the 
energy loss is a normalized distribution, n a defined in 
Eq. (21), which accounts for the theoretical expectation 
(-Btt, known as a Bichsel parameterization) and the reso- 
lution of the TPC for pions (ov). 



r V = log((dE/dx)/(B^))/a 7l 



(2) 



The resultant distribution in each transverse momentum 
range is fit with a six-Gaussian function (one per particle- 
species/chargc). The Gaussian widths are considered to 
be the same, independent of particle type, and single- 
Gaussian centroids are defined by the theoretical expec- 
tations constrained by the identified proton and pion 
measurements from topologically reconstructed weakly 
decaying particle yields 17j. Further details of thepar- 
ticle identification technique can be found in Refs. [l|, 0] • 
Raw data yields are corrected for single-track ineffi- 
ciencies evaluated via Monte-Carlo tracks embedded into 
real data events. We define single-track efficiency as 
the fraction of Monte Carlo tracks embedded into real 
Cu+Cu events that have been reconstructed. The ef- 
ficiencies are derived for each different event multiplic- 
ity bin and particle species. For high-px tracks (pt> 
2 GeV/c) in 200 GeV Cu+Cu events, the efficiency is 
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found to be 85% on average, with a weak (<10%) cen- 
trality and Pt dependence. In the analysis, pion and 
(anti)proton abundances are extracted from the n a dis- 
tribution using the finely calibrated centroid positions. 
The derived kaon yields are then smoothed to reduce sta- 
tistical fluctuations, using a Levy [l| fit. These assumed 
kaon spectra are then used for a final fit to determine 
the pion and proton yields. The systematic uncertain- 
ties from this procedure, on the spectra, are 2-10% for 
pions and 5-11% for protons, decreasing smoothly with 
Pt in the measured range. An analysis solving simulta- 
neous equations to assumed pion, kaon and (anti)proton 
distributions (bin counting), derived results that are 5- 
10% (5-20%) higher for pions (protons). This difference 
is the dominant systematic uncertainty on particle spec- 
tra. An additional systematic error of 5% resultant from 
the uncertainty in the single-particle efficiency determi- 
nation is added in quadrature. The total systematic er- 
ror for pion spectra ranges from 9% (at 3 GeV/c) to 13% 
(at 10 GeV/c). For protons, the error ranges from 21% 
(at 3 GcV/c) to 23% (at 6 GeV/c). These uncertainties 
are similar to the earlier evaluation from Au+Au data 
analysis Q. Systematic uncertainties from other possi- 
ble sources such as the momentum resolution (studied by 
embedding) and the uncertainty in determination of the 
centroid position (within the particle identification pro- 
cedure) are negligible. The corrected transverse momen- 
tum spectra for ir^ and (anti)protons at y / ijv r /v=200 GeV 
are shown in Fig. [T] (a) and (b) respectively The reach 
in transverse momentum is limited only by the avail- 
able statistics. Additional pr reach for pion identifica- 
tion is due to a larger separation from the kaon peak 
(Aa vK ~ 2cr, Aa pK ~ la). 

Figure [2] shows tt~ /tt + and p/p ratios in Cu+Cu data 
at y/SNN=200 GeV for the four centrality bins. The data 
show no systematic trends versus centrality within uncer- 
tainties, and a weak (if any) decreasing p/p with trans- 
verse momentum (as observed in Au+Au collisions 0]). 
Thus, to improve the statistical uncertainties in the fol- 
lowing discussion, data are averaged over particle charge. 

The spectral data alone can convey only a limited mes- 
sage. To delve into properties of the resultant data, ra- 
tios are taken. The first such ratio is termed the nuclear 
modification factor (Raa), defined in Eq. (TT]). We find 
that the pion spectra are suppressed in the most central 
(head-on) Cu+Cu data at Jj^=200 GeV (Fig.©. For 
the peripheral (glancing) collisions, a small enhancement 
is observed. To expose the features of the modifications 
of the hadron spectra in Cu+Cu and Au+Au collisions 
we study Raa as a function of the amount of matter 
participating in the collisions. For both systems Raa is 
evaluated within several fractional cross-section bins and 
as a function of the number of participating nuclcons. 
Figure [3] (a) shows the results of this comparative anal- 
ysis using the most central 0-12% (open squares) and 
mid-peripheral 40-60% (open circles) Au+Au data. For 
the most central events the suppression level is found to 
be different between the systems. The resultant spec- 



tra from Au+Au collisions are more suppressed than in 
Cu+Cu data. According to the Glauber calculation the 
mid-central (20-40%) Cu+Cu collisions (closed circles) 
and mid-peripheral (40-80%) Au+Au data (open circles) 
have similar numbers of participating nucleons (see Ap- 
pendix A for details). For this selection of ccntralities 
within the two systems we find that numerical values 
of Raa agree within the uncertainties. This agreement 
suggests a correlation of the suppression with the initial 
volume of the collision system. 

In Fig. [3] (6) we present the pt averaged Raa for pions 
(5<pt<8 GeV/c) as a function of the number of partici- 
pating nucleons calculated for both Cu+Cu and Au+Au 
collisions. The agreement between Au+Au (open circles) 
and Cu+Cu (closed circles) is striking and demonstrates 
that the nuclear modification factor for pions is a smooth 
function of the number of participating nuclcons (inde- 
pendent of the collision system) . 

Similarly, we explore the systematics of baryon produc- 
tion in Cu+Cu and Au+Au systems by comparing the 
Raa for protons and anti-protons. Figure [31(a) shows the 
Raa distributions averaged over p and p for four cen- 
trality bins of Cu+Cu events. The data at hand does 
not differentiate if collision volume (iV par t) or fractional 
cross-section effects are driving the high-py suppression 
for baryons due to the larger systematic uncertainties 
for (anti)proton measurements. Nevertheless, we observe 
that proton production in the peripheral Cu+Cu events 
is consistent with binary scaling expectations, and the 
suppression is setting in as one progresses from the pe- 
ripheral to central events. An overall similar central- 
ity dependence was observed between the Au+Au and 
Cu+Cu data at the same energy (see Fig. 2] (&)), albeit 
Cu+Cu integrated Raa values seem lower than the re- 
spective Au+Au data points. We emphasize that the sys- 
tematic errors are uncorrelated between the systems, and 
both measurements are similar within the experimental 
uncertainties. 

The similarity between the different systems at the 
same number of participants is also evident in other as- 
pects of the data at lower pr [H|. The smooth depen- 
dence of the nuclear modification factor could be inter- 
preted as a consequence of medium induced energy loss 
of partons traversing the hot and dense medium. For 
the smaller systems sizes, either peripheral Au+Au or 
Cu+Cu data, the path length traversed is smaller (on 
average) than for the larger system (central Au+Au). 
As observed in the data, a smaller energy loss is thus 
predicted @. 

Another dramatic effect observed in Au+Au data is 
the relative enhancement of protons to pions in the 
intermediate-pT region as compared to pp and e + + e~ 
collisions Q as well as for other baryon to meson ra- 
tios [l9j|. This enhancement is found to be strongly de- 
pendent on the centrality of the collision, as illustrated 
in Fig. \5\ The most peripheral A+ A data is shown to ex- 
hibit little or no enhancement in this ratio, with respect 
to pp collisions at the same energy. A similar increasing 
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FIG. 3: (Color online) (a) Nuclear modification factor, Raa, for charged pions (tt + + ir~) in Cu+Cu (filled symbols) and 
Au+Au (open symbols) collisions at y / ijvjv=200 GeV. Error bands are shown for most peripheral and most central Cu+Cu 
data to represent evolution of the systematic uncertainties for this dataset. Error boxes at Raa=1 represent Cu+Cu scale 
uncertainties due to the number of collisions and from pp spectra normalization, (b) Integrated pion Raa over the range 
5<pt<8 GeV/c versus iV par t. The bands represent the systematic uncertainty on ratios. An additional scale error due to pp 
normalization (~14%) is not shown. 



trend of favorable baryon production with centrality is 
observed in the Cu+Cu collision system. The peak of 
the enhancement is observed in the region py~2 GeV/c 
in Au+Au, at a slightly lower transverse momentum than 
the range measured in this analysis. At higher transverse 
momenta the enhancement over pp collisions diminishes 



to the level expected from vacuum fragmentation. 

The baryon to meson ratio (p+p) / in Cu+Cu 
and Au+Au collisions shows similar trends for an equiv- 
alent number of participating nuclcons. To further quan- 
tify this observation Fig. [5] (b) shows the proton to pion 
ratio (for hadrons with 3<£>t<4 GeV/c) measured in 
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FIG. 4: (Color online) (a) Nuclear modification factor, Raa, for protons and anti-protons (p + p) in Cu+Cu (filled symbols) 
and Au+Au (open symbols) collisions at y / ijvjv=200 GeV. Error band is shown for most central Cu+Cu data to represent 
characteristic systematic uncertainties for Cu+Cu data. Error boxes at Raa = 1 represent Cu+Cu scale uncertainties due to 
the number of collisions and from pp spectra normalization, (b) Integrated (anti)proton Raa over the range 5<pt<6 GeV/c 
versus N pal t- The bands represent the systematic uncertainty on ratios. An additional scale error due to pp normalization 
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FIG. 5: (Color online) (a) Ratio of protons and anti-protons to charged pions, versus transverse momentum, for Cu+Cu 
collisions at y / ijvjv=200 GeV. For clarity, only one systematic error band is shown for Cu+Cu data (most central events), 
uncertainties on data in other centrality bins are similar in magnitude. Systematic errors for Au+Au data are not shown, for 
details see 0]. (b) Average (p + p)/(ir + +iv~) ratio for 3<pt<4 GeV/c is shown as function of centrality (iVp ar t) for Cu+Cu 
and Au+Au data. Error band represents the pp measurement with uncertainty. 



Cu+Cu and Au+Au collisions as a function of iV part . We 
find that this ratio is also sensitive to the initial volume 
of the collision system and exhibits the same quantitative 
IVpart dependence irrespective of the collision system. 

As discussed earlier, it is found that in the kinematic 
range of our measurements baryons are produced pre- 



dominantly from gluon fragmentation [20(. It is thus 
expected that an increase in the baryon to meson ra- 
tio in the intermediate- to high-px range would be re- 
lated to gluon sources. To explain the presented data 
one could consider, for example, that a gluon jet could 
be more easily propagated through the medium than a 
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quark jet, leading to an increase in the number of pro- 
tons in the intermediate-py region. This, however con- 
tradicts theoretical predictions where an opposite effect 
was expected Alternatively, more gluon jets could 
be initially produced, or induced (for example, in the ra- 
diative energy loss scenario), for the more central data. 
The latter appears to be the more plausible, as the high- 
est pt data exhibits little or no enhancement over the 
pp data, indicating a similar energy loss for gluons and 
quarks (see Fig. [5]). Alternative approaches to explain 
the phenomenon observed in the data, have also been de- 
veloped. For example, the recombination/fragmentation 
picture of thermal/shower partons has had success at de- 
scribing this in Au+Au data [3]. Further information 
on the relative energy loss of quark and gluon jets can be 
extracted from the data by comparing the nuclear modifi- 
cation factors of proton and pion data (Figs. [3] and [4]). At 
high-pT (above 5 GeV/c), however, the two suppression 
factors are found to be the same within the systematic 
uncertainties, suggesting a similar energy loss of quark 
and gluon jets in Cu+Cu collisions. 

In conclusion, new results on high-pT identified pion 
and proton spectra are presented for several centrality 



bins in Cu+Cu collisions at ^/sjvat=200 GeV. The data 
are found to exhibit similar systematic trends over a wide 
range of transverse momenta as Au+Au collisions at the 
same energy with a similar number of participants. The 
suppression pattern observed versus the number of par- 
ticipants in Au+Au data is followed by the Cu+Cu data 
to a large degree. The participant coverage in these 
Cu+Cu collisions is in a region where the suppression 
effects are turning on. A detailed study of the proton 
to pion ratio reveals similar systematic dependencies to 
that found in Au+Au data. Specifically, the increase in 
proton yield at intermediate transverse momenta persists 
for the much smaller Cu+Cu system. 

Further studies have shown similar suppression of pro- 
tons and pions at high-pT- Within the context of the 
connection between the detected pions and protons and 
quark and gluon jets suggested in the introduction, these 
results indicate similar partonic energy loss for both glu- 
ons and quarks. The amount of energy loss suffered by 
the partons is found to be -/V par t dependent. Within the 
experimental uncertainties, the suppression for different 
collision species is found to be invariant for the same 
number of participants. 
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APPENDIX A: MONTE CARLO GLAUBER 
MODEL RESULTS FOR THE CENTRALITY 
BINS USED IN THE PAPER 



TABLE III: Number of participants A par t and number of 
binary collisions Abin from the Monte Carlo Glauber model 
calculations for 200 GeV central triggered Au+Au collisions. 



TABLE I: Number of participants Apart and number of bi- 
nary collisions Abin from the Monte Carlo Glauber model cal- 
culations for different centrality bins of minimum bias Cu+Cu 
collisions at 200 GeV. 



Centrality bin A p 



A bin 



0-12% 



315.7+^j 900.3lgg j 



Centrality bin A part 



A bin 
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74.6t 1 1 ;2 123.6t»;$ 



1+4.2 



45.9tT 6 62.9± 3 . 7 



21.5] 



22.7+j-i 



TABLE II: Number of participants A part and number of bi- 
nary collisions Abin from the Monte Carlo Glauber model cal- 
culations for different centrality bins of minimum bias Au+Au 
collisions at 200 GeV. 
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